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A search for double beta decays of tin isotopes with enhanced sensitivity 
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A search for the various double beta decay modes of ^^''Sn and ^^^Sn has been performed on 
75kg-days of data. New half-life limits for excited states in ^^^Sn have been obtained including a 



lower limit for decay into the first excited 2+ state of ^^"Te of Ti/a > 0.87 x 10^" yrs (90%CL) 
and into the first excited O"*" state of T-1/2 > 1.08 x 10^" yrs {90%CL). Ground state and ex- 
cited state transitions of '^^'^Sn have also been experimentally explored. A limit for the 2i'EC /EC 
of Ti/2 > 1.8 X 10 yrs (90%CI/) is obtained. The non-observation of de-excitation gammas 
from the 0+ at 1888.5 keV results in a lower half-life limit on the OuEC/EC decay of "^Sn of 
2^1/2 > 0.8 X 10^^ yrs (90%CI/), despite a possible resonant enhancement of the decay rate due to 
degenerated states. 

PACS numbers: 



I. INTRODUCTION 

The Standard Model of Particle Physics has been ex- 
tremely successful but is not thought to be the final the- 
ory. The quest for physics beyond the Standard Model is 
gathering pace, with searches performed at accelerators 
such as Tevatron, HERA and soon at the LHC. There 
are also non-accelerator experiments, such as the search 
for the rare process known as Neutrinoless Double Beta 
Decay. 

The classical Double Beta Decay is allowed by the 
Standard Model. In this reaction two neutrons in the 
same nucleus decay simultaneously to emit two anti- 
neutrinos and two electrons. This possible reaction was 
first discussed by M. Goeppert-Mayer[l| in the form of 

{Z,A) ^{Z + 2,A) + 2e- +2i7e (2i//3/3 - decay). (1) 

As this is a higher order process, practically to observe 
this requires an isotope in which single beta decay is 
forbidden or strongly suppressed. The energy spectrum 
from the two electrons has similar features to that from 
single beta decay, a continous spectrum ending at a well 
defined end-point determined from the Q- value of the re- 
action. 

The classic papers of Majorana [3], in which he dis- 
cussed a two-component neutrino, led Racah [3| and 
Furry to first discuss a neutrinoless mode in the form 
of i 

{Z,A) ^{Z + 2,A) + 2e- (Oj^/3/3 - decay). (2) 

The signature of this mode is very different to the clas- 
sical case. Here the summed energy of the emitted two 
electrons is exactly equal to the Q-value of the reaction. 
The experimental concept therefore is to observe an iso- 
tope capable of decaying via double beta decay and to 
search for a peak at the endpoint of the spectrum. 



This process clearly violates total lepton number con- 
servation by two units and so is forbidden by the Stan- 
dard Model. This is in contrast to neutrino oscillations 
which violate individual flavour lepton number but con- 
serve total lepton number. For this decay to occur neu- 
trinos must be massive Majorana particles, however any 
AL = 2 process can contribute to the decay. The most in- 
teresting mode is the light Majorana neutrino exchange. 
Here the measured quantity is known as the effective Ma- 
jorana neutrino mass which is given by 
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which for the case of CP-invariance, e^*"' = 0, tt, is 

(m,J =1 TOit/Ji ± maC/la ± "^af/es I (4) 

For double beta decay the decay rate and thus inverse 
half-life is strongly dependent on the available Q-value. 
In the classic case, 2i/f3/3-deca.y, the rate is proportional to 
Q^^ whilst the OvP/S - decay mode scales with Q^. For 
this reason, only high Q-value decays (typically above 
2 MeV) are currently worth studying. This criterion re- 
stricts the 35 candidate isotopes to 11. 

There are also double beta decay modes (both 
21^(313 - decay and Oz^/3/3-decay) in which gamma rays are 
emitted in addition to the two electrons, thus giving 
these modes distinct signatures. In principle, a finely seg- 
mented detector could observe the two electrons in one 
segment and simultaneously the emitted gamma ray in 
another segment. Since the energies of the electrons and 
gamma ray are known this gives a strong experimental 
signature which is distinct from the background. Study- 
ing these modes may shed some light on new physics 
and provide some new information for nuclear matrix el- 
ement calculations. Despite their reduced Q-values there 
are benefits in studying these modes. 



These modes come from two sources, double beta de- 
cays to excited 0+ and 2+-states of the daughter and 
double positron decays in which positrons are emitted 
instead of electrons. Observations of 2i//3/3 - decay tran- 
sitions into the first excited 0+-state for ^'^°Mo and ^^'^Nd 
have already been reported @, H S Q • 

For the double positron decays, combinations of elec- 
tron capture (EC) and positron emission can occur: 

{Z,A)^{Z-2,A) + 2e+ + {2y,) (3+13+ (5) 

e- + {Z,A)^{Z-2,A) + e+ + {2ve) 13+ 1 EC (6) 

2e- + {Z,A)^{Z -2,A) + {2v,) EC / EC (7) 

In particular, the (3+ / EC mode shows an enhanced sensi- 
tivity to right handed weak currents [9|. For each positron 
there is the possibility of also observing one or both of the 
511 keV photons. These modes therefore can provide an 
extremely clean signature of up to five coincident energy 
deposits. However, for each generated positron the avail- 
able Q- value is reduced by 2meC^, which leads to much 
smaller decay rates than in comparable 0vP(3 - decay. 
For (3+(3+-deca.y to occur, the Q-values must be at least 
2048 keV. Only six candidate isotopes are known to have 
such a high Q- value. 

The full Q- value is available only in the EC/ EC mode 
but is difficult to detect. In a neutrinoless EC /EC to the 
ground state of the daughter, a monoenergetic internal 
bremsstrahlung must be emitted and requires electron 
captures of both K- and L-shells. If the initial and final 
states are degenerate in the context of radiative EC /EC 
then a resonant enhancement in the decay rate could oc- 
cur 10]. The signal observed would come from the de- 
excitation photons. 

The three combinations of electron capture and 
positron emission have not yet been observed, not even in 
the neutrino accompanied mode, though there has been 
a weak indication of the observation of ^'^"Ba decay in 
geochemical experiments [ll| . 

There is a wide variety of possible gamma ray emis- 
sions. This paper explores those double beta transitions 
of tin isotopes which emit gamma rays. There are three 
double beta isotopes of tin, ^^^Sn and ^^''Sn in the two 
electron mode and "^Sn for (3+ /EC and EC / EC de- 
cays. The Q-values of the transition for each of the three 
isotopes are 366 keV, 2287 keV and 1922 keV and the 
natural abundances are 4.63%, 5.79% and 0.97% respec- 
tively. As there is no excited state of interest for ^^^Sn 
decay we focus on the decays 



i24Sn 
2e- + "^sn 



i24Te -F2e- + (2Pe)+7 (8) 
"2Cd +(2z/,)+7 (9) 

"2Cd +e+ + (2i/e)+7 (10) 



^^"^Sn is one of the eleven isotopes with Q-values larger 
than 2MeV yet there is at present no proposal for a 
large scale experiment and few calculations. There are 



no theoretical predictions of half-lives for -'^^^Sn. There 
are calculations reported for 2^(3(3 - decay ground state 
transitions in ^^''Sn, varying from 0.7 — 3 x lO^^yrs and 
estimates for the first excited states 0+ and 2"'"-states are 
2.7 X 10^1 yrs and 6 x lO^^yrs respectively. [ll,[il 

Calculations within the single state dominance model 
for the 2v 13+ / EC and 2v EC /EC ground state tran- 
sitions in ^-'^^Sn have been made|19i|. The half-lives ex- 
pected are of the order of 10^^ yrs {EC /EC) and 10^'* yrs 
{(3+ /EC). However it is possible that an enhancement 
of the rate for the (3+ /EC modes occurs due to right- 
handed weak currents. 

The radiative neutrinoless EC /EC rate to the ground 
state of ^^^Cd has been estimated to be of the order of 
10~^^yrs~"'^. This rate may be enhanced by a factor of 
10^ by a possible resonance condition [1(||. The char- 
acteristics of this mode are the monoenergetic internal 
bremsstrahlung gamma at Q — Ek — E^ = 1888.5 keV, 
with Ek,El the K- and L-shell binding energy of ^^^Cd. 
There could also be a degeneracy between ^^^Sn and an 
excited 0+ state at 1870.9 keV in "^cd, which would ful- 
fil the resonance enhancement condition for neutrinoless 
EC /EC 0]. 

Two recent studies on tin decays produced half-life lim- 
its for various decay modes in the region of 10^^ — 10^^ yrs 
[20l . |2l| . The aim of this measurement was to improve on 
those limits, if possible, by another order of magnitude. 
For this analysis, gamma line energies were taken from 
[23| and the latest atomic mass determinations from [2J] . 



II. EXPERIMENTAL SETUP 

The gamma-spectrometric measurements were per- 
formed in the Felsenkeller underground laboratory near 
Dresden (Germany) at a depth of 110 mwe. 36 tin bars 
of 10.6 cm^ area and 4 mm thickness plus 4 tin sheets of 
24.2 cm^ X 3 mm were placed in a Marinelli beaker so that 
the endcap of the measuring detector was surrounded by 
the material. The total mass of tin used was 1.24 kg. 
A passive shielding consisting of 170 mm lead of graded 
qualities protected from the ambient background radia- 
tion in the measuring chamber. ^^^Rn concentration in- 
side the shielding was reduced by flushing with gaseous 
nitrogen. The gamma spectrometer used a well-type low 
background p-type HPGe detector with 150 cm^ sensi- 
tive volume (30 % relative efficiency, FWHM at 1.3 MeV: 
2.0 keV). For further background suppression a U-shaped 
cryostat configuration was applied. Since the muon flux 
in the laboratory is reduced by a factor of 50 compared to 
the surface value, no additional active shielding was nec- 
essary. The pulse processing electronics used a standard 
spectrometric amplifier and the spectrum acquisition was 
performed by an 8 k channel MCA. Spectra of the tin 
sample were accumulated during two measuring periods 
of 598 and 853 hours duration, respectively. Thus, the 
total exposure amounts to 74.96 kg days (where the mass 
refers to natural tin). A blank spectrum of the empty 



shielding was measured for 762 hours. 
Just before data collection commenced, the Germanium 
detector was exposed to a mixed source of ^"^^Am, ^^^Cs 
and ^°Co. The gamma lines observed from this source 
were fitted with gaussian functions. The fitted means 
were used to accurately calibrate the ADC channel to 
energy relation and the widths were used to derive a func- 
tion for the energy resolution (values given in Tab.|T|. 
The efRciency and response of the setup was studied us- 
ing detailed GEANT4 based Monte Carlo simulations, 
which modelled all parts of the germanium crystal, the 
crystal housing, the tin pieces and the background shield- 
ing. Samples of 10^ gammas were simulated inside the 
tin volume for each energy of interest and the number 
of events in the photo-peak recorded in the crystal was 
used to determine the efficiency for observation, along 
with Poissonian uncertainties. Furthermore, know con- 
taminants were simulated to confirm understanding of 
the shape of the background distribution. 
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FIG. 1: Total spectrum (black), background spectrum (red). 
A clear bremsstrahlungs contribution at lower energies is seen. 



III. DATA 



The obtained spectrum compared to a background run 
is shown in Fig. [TJ Due to the better low energy be- 
haviour of the detector in the Felsenkeller with respect 
to the one used in [21] the 46.5 keV line from ^-'^"Pb be- 
comes visible, showing the presence of lead in the tin. 
2iOpb decays with a half life of 22.20 years via ^^^Bi 
(5.012 days) and ^lOpo (138.38 days) to stable ^o^Pb. 
Simulations were used to verify that the beta emitter 
^^°Bi {E,nax = 1162.1 keV) is responsible for the brem- 
strahlungs background below 600 keV. ^-'^"Po undergoes 
a-decay but emits a weak gamma line at 803.10 keV 
which is visible in the tin spectrum. Further contami- 
nation of ^^'*Pb and ^^''Bi (^^^Rn successors) as well as 

of 228Ra, 228rpi^^ 40j^^ IST^g ^^^^ QOq^ ^^^ ^jgjble in the 

tin spectrum at count rates below 0.2 counts/hr (Fig[2]). 

Weak saw-tooth like features are observed in both tin 
and background spectra at 693 keV and 596 keV. These 
can be attributed to '''^•''^Gc{n,n') interactions [22] ■ The 
same features were observed more strongly in data col- 
lected in a surface laboratory, where the neutron back- 
ground is expected to be higher, for Ref. [21]. These pre- 
vious data were used to determine the functional form: 
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FIG. 2: Background subtracted and thus the spectrum of the 
tin sample. 



IV. RESULTS 

The search performed relies on decays of "^"^^Sn and 
^^^Sn. The relevant experimental data, including energy 
resolution, efficiency and the background at the energy 
of each line searched for are compiled in Tab HI 



A. The ^^*Sn system 



f{E) = K.(1.0-0.04(S-So)) 

when (E> Eq) & (1.0 
f{E) = Otherwise 



(11) 
0.04(£; - So)) > 



which was found to describe the shape of these features 
well. Here, So is the starting point of the feature, taken 
to be 692 keV and 595 keV respectively, whilst k is the 
magnitude parameter, which varies with the neutron flux. 



A Bayesian analysis approach 25], similar to that de- 
scribed in Ref. [21] was adopted. The most likely values 
for the magnitude of a linear background component and 
the amplitude of a gaussian peak at fixed gamma energy 
were determined from a binned maximum likelihood fit. 
The mean and width of the gamma-line under investiga- 
tion were fixed in the fit, which was carried out over a 
range of ±30cr around the known gamma energy. When 
gamma lines from known radioactive contaminants fell 
within this fit range, their amplitude was also included as 



Energy (keV) 


AE (keV) 


e(%) 


B (counts/keV) 








Q 


/? 


553.8 


1.739 


3.51 


476 


-0.62 


602.7 


1.763 


3.36 


425 


-0.49 


606.5 


1.765 


3.38 


477 


-0.60 


617.3 


1.770 


3.34 


533 


-0.69 


694.7 


1.805 


3.06 


295 


-0.29 


695.0 


1.806 


3.06 


303 


-0.30 


713.8 


1.814 


3.01 


418 


-0.48 


722.9 


1.818 


3.00 


377 


-0.42 


815.1 


1.858 


2.76 


167 


-0.13 


851.1 


1.873 


2.67 


139 


-0.10 


1054.0 


1.953 


2.31 


118 


-0.07 


1253.4 


2.022 


2.05 


126 


-0.07 


1488.9 


2.093 


1.82 


166 


-0.10 


1888.5 


2.184 


1.49 


123 


-0.06 



TABLE I: The energy resolution, AE (FWHM), gamma de- 
tection efficiency, e, and fitted background, B, for each peak 
energy. The magnitude of uncertainty in e ranges from 0.01- 
0.02%. The background continuum is defined by the poly- 
nomial B — a + PE and gives the total counts over the 
1451 hours of data collected. 
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FIG. 3: Level diagram of Sn. The energy of the dominant 
mode decay gamma is given for each level along with the 
percentage branching ratio in brackets. 



a fit parameter (again energy anA resolution were fixe(i). 
Thie effect of the ^^'^^Ge(n,n') interactions at 693 keV and 
596 keV was also included, when relevant, by adding the 
function in Eqn. Ill[ allowing the parameter n to vary as 
a parameter of the fit. 

The fitted signal amplitude and uncertainty deter- 
mined from the likelihood fit, 9g ± Sg was then used to 
derive a 90% lower limit on the half-life for each decay 
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FIG. 4: Peak range around 602.7 keV. Natural radioactivity 
from ^°®T1 and ^^''Bi produces the lines at 583 and 609 keV 
respectively. 



mode using Eqn. 1121 
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half 
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Thali > 



{0s 



1.28(5,) 
iiivce'01n2 



(Os > 0), 



ISO ''live 



i0s < 0) 



(12) 



(1.28(5,) 

Here A'iso is the number of candidate nuclei in the tin 
for the given decay, tuvc is the duration of data collec- 
tion in years, e is the efficiency for observing the given 
gamma signal determined from simulations and ip is the 
branching ratio for the gamma line in question. The fac- 
tor of 1.28 gives one-sided 90% limits. For each fit, the x^ 
goodness of fit was determined and in all cases indicated 
good agreement with the data. 

The level scheme of -'^^''Te (given in Fig. [3]) shows that 
the higher level decays all dominantly proceed via the 
intermediate 2f-state at 602.7 keV. Hence, studies of a 
possible gamma peak at 602.7 keV provide information 
on all these decays. The region of interest about this 
energy is shown in Fig. |4] along with the fitted signal 
gaussian and backgrounds. The efficiency for observing 
gammas of this energy from the tin was determined to be 
3.36±0.02% from simulations and the calculated energy 
resolution at this point is 1.76 keV (FWHM). With 1.4 ± 
4.6 possible events in the gaussian peak a half-life of 



y^^^0(2).(Q4 



2+(602.7fcey)) > 8.7 (11.0) x lO^Vs 
(90 (68) %CL) (13) 



was obtained, an improvement of a factor 28 on existing 
limits [2l| which can be attributed to the larger size of 
germanium detector and the reduced background due to 
the underground location of this experiment. 

Although decays into higher excited states will pro- 
duce two or more gammas, it is unlikely that they will 
all be emitted in the direction of the detector, and there- 
fore, studies of higher states were performed by search- 
ing separately for the accompanying gammas. The fitted 
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FIG. 5: Level diagram of ^^^Sn. The energy of the dominant 
mode decay gamma is given for each level along with the 
percentage branching ratio in brackets. 



events and (derived half-life limits for these decays are 
compiled in Tab. [TTl The branching ratio for the de- 
excitation chain, which is given in Fig. [3] and can be less 
than 100% for the higher excited states, was included in 
the calculation of half-life. In some cases the best limit 
on half-life is derived from one of the lower de-excitation 
steps in the decay chain (indicated by a f in the table). 



Excited state 


Gamma energy (ies) 


Events 


Tij-i 




(keV) 




(10^° yrs) 


2+ (602.7) 


602.7 


0.6 ± 12.3 


0.87 


0+ (1156.5) 


553.8, 602.7 


-44.4 ± 10.6 


1.08 


2+ (1325.5) 


722.9, 602.7 


6.5± 8.8 


0.62t 


Oj (1656.7) 


1054.0, 602.7 


-1.8 ± 6.7 


1.13 


Oj (2020.0) 


695.0, 722.9, 602.7 


18.1 ±12.1 


0.38t 


2+ (2039.3) 


713.8, 722.9, 602.7 


-8.1 ± 8.5 


0.62t 


2+ (2091.6) 


1488.9, 602.7 


-9.0 ± 4.7 


1.08 



TABLE II: Half-life limits (90 % CL) for all possible excited 
0^,2^-state transitions (Oi^ and 2v) of ^^^Sn. Given are the 
excited states, all possible gamma lines with the dominant 
mode searched for in bold, the total number of fitted events 
in the peak and deduced half-life limits. For decays marked 
I the better half-life limit of the 602.7 keV gamma in the first 
row applies. 



B. The "^Sn system 

The second system to analyse is ^"'^^Sn , which has the 
decay level scheme shown in Fig. [5] It offers a more 
complex search pattern as in addition to the de-excitation 
gammas more photons can be emitted in the EC / EC or 
(3^ /EC process. 
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FIG. 6: Data in the range around 617.3 keV. The line at 
609.3 keV stems from ^^'*Bi decay in the natural ^^*U decay 
chain. 



Again, all decays into higher excited states proceed 
via a 2^ state, this time at 617.3 keV. Thus, a search 
for this line places a half-life limit on all excited state 
transitions, taking into account the branching ratios of 
the de-excitation. Fig. [6] shows the fitted region where 
no peak was observed, resulting in a half-life limit of 



rp 0(2)uEC/EC+[3+/EC /r,+ 
-'1/2 I" 



2+(617.6 keV)) > 
1.8 X 10^^ yrs (90%CL). (14) 



Decays to excited states were treated in a similar manner 
to the ^^^Sn system. The results of searches for all the 
potential gamma lines and the derived half-life limits are 
compiled in Tab. IIIII 

Decay to the 0~''-state at 1870.9 keV is particularly 
interesting as this state could be degenerate with the 
ground state of ^^^Sn , which could result in a resonantly 
enhanced neutrinoless EC/EC rate. De-excitation of 
this state is dominated by the emission of a 1253.4 keV 
and 617.3 keV gamma. No obvious peak was observed at 
either of these energies, resulting in a half-life limit on 
this decay mode of 



ji 0(2)^'BC/BC('Q^ 



0+(1870.9 keV)) 



-|19, 



> 1.8 X lO'-yrs 
{90%CL). (15) 



The P^ /EC modes can only populate excited states 
up to Q-1.022MeV, which corresponds to levels below 
900 keV for "^Sn. Only the 2+-state at 617.3 keV sat- 
isfies this criteria. Therefore the half-life limit from the 
617.6 keV hue search for EC /EC, given in Eqn. [T4l also 
applies for the Q{2)v (3^ /EC -decay modes into the first 
excited 2^ state. 

A search was also performed for the Qv EC /EC ground 
state transition, based on the emission of an internal 
bremsstrahlung photon [2q . Due to energy and momen- 
tum conservation this requires one K- and one L-shell 



Excited state 


Gamma energy 


Events 


T\/2 




(keV) 




(10'^ yrs) 


2t (617.3) 


617.6 


-5.4 ± 10.1 


1.8 


0+ (1224.1) 


606.5, 617.6 


-1.1 ± 10.4 


1.8 


2+ (1312.3) 


694.7, 617.6 


15.2 ± 12.1 


0.5t 


0+ (1433.2) 


815.8, 617.6 


8.3 ± 8.0 


0.7t 


2+ (1468.7) 


851.1, 617.6 


-0.5 ± 7.4 


1.3t 


O;^ (1870.9) 


1253.4, 617.6 


-5.0 ± 5.8 


1.7t 



TABLE III: Half-life limits (90 % CL) for excited state tran- 
sitions (0:^ and 2v) of ^^^Sn. Given are the excited states, all 
possible gamma lines with the dominant mode searched for 
in bold, the total number of fitted events in the peak and de- 
duced half-life limits. For decays marked f the better half-life 
limit of the 617.3 keV gamma in the first row applies. 



which resulted in a half-life Hmit of 



Ti/2 > 2.0 X lO'^^yrs (90%CL) 



(16) 



due to the non-observation of a peak. 
Uncertainties in the atomic masses mean that the peak 
position is only know to about ±2keV [2jy. Therefore, 
the peak position was systematically varied by 0.5 keV 
between 1886 and 1890 keV and the fit was repeated for 
the corresponding position. The conservative lower hmit 
for this decay is taken as 



rp OvEC/EC i(^+ 
J 1/2 (0 



0+ ) > 0.8 X lO'Vrs (90%CL) 



the worst half-hfe hmit found in this range. 

V. SUMMARY AND CONCLUSIONS 



(17) 
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FIG. 7: Data around 1888.5 keV, and the fitted peak limits. 



Double beta decay, in its various forms, is a power- 
ful tool for investigating neutrino properties. Although 
it has received little interest in the field so far, tin 
presents two isotopes of interest: ^^^Sn which is suitable 
for Qv(3~ I3~ searches and ^^^Sn available for decay via 
the (3^ /EC mode. In this paper, new half-life limits are 
presented for various decays of these two isotopes, based 
on the non-observation of characteristic gamma- lines. All 
of the limits presented here improve on previous values 
given in Ref. 21| by an order of magnitude or more. With 
limits of order 10^" yrs for a number of decay modes, this 
paper might trigger some new theoretical efforts to cal- 
culate 2vPj3 - decay half-lives in the tin systems. During 
submission of this paper we became aware of a similar 
study in Ref. [13|. 



capture. Additionally the bremsstrahlung photon has to 
be monoenergetic, with an energy of the Q- value reduced 
by the K- and L-shell binding energy of the daughter. In 
the case of ^"'^^Sn this implies a gamma of 1888.5 keV. The 
calculated detection efficiency for such a photon is 1.491 
± 0.013 %. Fig.[7]shows the energy region for this search. 
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